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ABSTRACT 


The primary objective of this effort consisted of a detailed study 
of the history of the motion of the moon. More precisely this grant has 
resulted in the support of David Rubincam in the "preparation of "The 
Early History of the Lunar Inclination". This has consisted of the 
completion of his thesis within a cooperative program between the 
University of Maryland and the Goddard Space Flight Center. In addition, 
several analyses were developed which are related to the determination 
of the effect of various refractive phenomena on the accuracy of measurements 
made through the earth's atmosphere. 
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I., GENERAL DISCUSSION 


The primary effort supported by this grant consisted of a research 
study, conducted by David Rub incam, entitled "The Early History of the Lunar 
Inclination^^, This study and the resultant thesis was in the form of a 
cooperative program directed by Dr, John 0^ Keefe of the Goddard Space Plight 
Center and the Principal Investigator at the University of Maryland, The 
objective of this study was to evaluate the lunar motion by numerically 
integrating lunar motion backward in time to an epoch T^en the moon was 
much closer to the earth and the inclination of the orbit was much greater 
than the current value. At this time the value and the role of viscosity 
of the earth in order to correct errors which had been made by investigators 
using earlier analyses. These errors had lead to improper conclusions on 
the mode of formation of the earth-moon system. 

The full description of this work and a discussion of the results 
is contained in the thesis presented to the University of Maryland by 
Mr. Rub incam. In addition, this material is also presented in a Goddard 
Space Plight Center Technical Document #X-592-73-328. 

Rather than repeat this discussion or attach the thesis in full to 
this final report, certain particularly relevant sections are presented in 
the next section. These include the abstract, the acknowledgments, and 
the table of contents which indicates the full range of topics addressed 
within the thesis. The introduction illustrates the 'basic method of attack. 
The discussion section contains the technical conclusions from the mathe- 
matical treatment and the implications on the formation of the earth-moon 


system. 
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THE EARLY HISTORY 
OF THE LUNAR INCLINATION 

David Parry Rubincam 



ABSTRACT 


The effect of tidal friction on the inclination of the lunar orbit to the earth's 
equator for earth-moon distances of less than 10 earth radii is examined. The 

results obtainea bear on a conclusion drawn by Gerstenkorn and others which 
has been raised as a fatal objection to the fission hypothesis of lunar origin, 
namely, that the present nonzero inclination of the moon's orbit to the ecliptic 
implies a steep inclination of the moon!^ orbit to the earth's equatorial plane in 
the early history of the earth-moon system. This conclusion is shown to be valid 
only for particular rheological models of the earth. In the case of a viscous 
earth, the results indicate that the problem of wrenching the moon out of an 
equatorial orbit into an inclined orbit to account for the present tilt of the lunar 
orbit to the ecliptic must be faced in the accretion theory of the moon's origin 
and possibly the capture theory, as well as in the fission theory. In this respect 

all three theories are on the same footing. A solution to the inclination problem 

1 

is presented. 

The treatment of tidal friction adopted here employs the approach of George 

Darwin and pursues his suggested solution to the inclination problem in great 

detail. The earth is assumed to behave like a highly viscous fluid in response to 

tides raised in it by the moon. The moon is assumed to be tideless and in a 

circular orbit about the earth. The equations of tidal friction are integrated 

numerically to give the inclination of the lunar orbit as a function of earth-moon 

lii 



distance. It Is found that if the radius of the lunar orbit is greater than 3.83 
earth radii, then the inclination of the moon's orbit to the earth's equator will 
increase if the moon is perturbed from an equatorial orbit, provided the earth's 
viscosity is greater than 10^® poises. The present inclination of the lunar orbit 
to the ecliptic can be explained if the moon's orbit is perturbed about 3° out of 
the equatorial plane at 3.83 earth radii, provided that the earth's viscosity is 
not less than 10^® poises. It is also found that if the viscosity is large (greater 
than 10^® poises), then, under certain conditions, the radius of the moon's orbit 
may actually decrease temporarily, and then increase; and further, that an 
upper limit can be placed on the inclination of the lunar orbit to the earth's 
equator when the moon is 3.83 earth radii distant from the earth, regardless of 
the moon's prior history. ^ 
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INTRODUCTION 


Men have speculated about the tides for centuries . An ancient Chinese 
scholar suggested that the earth lived, that the ocean was its blood, and that the 
tides were the beating of the earth's pulse (Darwin 1962, pg. 76). An Arabian 
scholar explained the rising of the tide as being caused by the heating of the 
ocean by sunlight and moonlight (Darwin 1962, pp. 77-79). Others also suggested 
that the tides were somehow caused by the sun and moon (Darwin 1962, pp. 79- 
85), but it remained for Isaac Newton to advance the correct explanation for the 
cause of the tides (Newton 1966). Newton realized that the lunar tides were 
caused by a combination of gravitational pull and centrifugal effects which would 

make the water in the oceans collect on the sides of the earth directly under and 

\ 

directly away from the moon, thus giving the earth a bulge. A similar argument 
holds for the solar tides. 

Newton's theory of the tides was carried forward, notably by Bernoulli, 
Laplace, Darwin, and Kelvin (Darwin 1962, pp. 86-88) to explain the rise and fall 
of the oceans on the earth. Their efforts culminated in the work of Doodson and 
Proudman (Doodson 1958). 

George Howard Darwin, son of the famous Charles Darwin, considered not 
only the problem of the tides raised on the earth by the moon, but also a more 
subtle problem: the action of the tides on the motion of the moon (Darwin 1880). 
He included in his investigations not only ocean tides, but also tides raised in the 
bulk of the earth as well; these latter tides are called earth or body tides. 

Darwin recognized that friction attending the tides, whether they are raised in 
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the oceans or in the earth, would have profound effects on the moon's orbit. In 
fact, tidal friction dominates the secular change in the moon's orbital elements. 

Darwin assumed the earth to be a homogeneous, incompressible viscous 
fluid in which tides were raised by the moon; the moon itself was assumed to be 
a point-mass. He expanded the tidal disturbing function in a Fourier series and 
integrated the equations for the secular change of the moon's orbital elements 
backwards in time in an attempt to uncover the past history of the moon. Darwin 
found that the moon orbited very close to the earth at some time in the distant 
past. He speculated that the earth and moon were once a single primitive body, 
and that resonance vibrations set up in the body by the sun caused the body to 
fission, thus throwing the moon into orbit about the earth. Tidal friction then 
caused the moon to move away from the ear^to its present distance. 

Jeffreys (1930) found that dissipation in the primitive body would be so 
great that the vibrations would be damped out, making it impossible for the 
moon to be torn out by the action of the sun. The fission theory of the origin of 
the moon fell out of favor. It has been reproposed in more recent times by 
Cameron (1963), Wise (1963), and O'Keefe (1969). 

Modern interest in tidal friction was rekindled by Gerstenkorn (Alfven 1963), 
who invoked tidal friction in a new hypothesis of lunar origin: capture of the 
moon. Gerstenkorn !s analysis lead him to propose that the moon was once an 
independent planet in an orbit which carried the moon close to the earth. The 
tidal interaction between the moon and earth captured the moon in a retrograde 
orbit, which subsequently flipped over into the prograde orbit we see today. 
MacDonald (1964) supported a many-moon hypothesis of lunar origin to 
overcome a time-scale difficulty in tidal evolution. Singer (1968) inves- 
tigated the problem of prograde capture, while the analysis of Goldreich (1966) 
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lead him to favor accretion of the moon from a swarm of particles in orbit about 
the earth. 

The problems associated with the intimately connected questions of the 
moon’s origin and its orbital evolution are seen to be of surpassing interest 
today. We will investigate here one of the problems of the early history of the 
moon: the inclination of the lunar orbit. 



CHAPTER V 


DISCUSSION 

A. Critique of Assumptions Made 

We shall now examine the important assumptions made in obtaining our 
results for strong tidal friction. 

One important assumption we have made is that the orbit of the moon re- 
mains circular throughout its history, i.e. the eccentricity e of the lunar orbit is 
zero. The work of Darwin (1880, Section V), Singer (1968), and MacDonald (1964) 
shows that weak tidal friction decreases the eccentricity as we look back into the 
past until the moon reaches about 3 earth radii from the earth, where the eccen- 
tricity undergoes rapid changes. Since the present value of the eccentricity is 
0.055, this would imply that neglect of the eccentricity when the moon was at the 
reference distance of 3.83 earth radii would lead to negligible error in con- 
sidering weak tidal friction. 

However, use of Darwin's treatment of the eccentricity for viscosities 
greater than 10^’ poises indicates that e increases with time until the moon 
reaches about 16 earth radii distance from the earth; at larger distances the 
eccentricity rapidly decreases. This indicates that the eccentricity could have 
been large for earth-moon distances of less than 16 earth radii. However, a 
nearly circular orbit for the moon over its whole history is by no means ex- 
cluded. The earth could have behaved as though it had a large viscosity when the 
moon was less than- 16 earth radii from it; beyond 16 earth radii the earth could 
have behaved as though it had a small viscosity. If this were the case, then if the 


57 



58 


moon were In a nearly circular orbit at 3.83 earth radii, the eccentricity would 
slowly grow to its present value as the moon moved outward to its present 
distance. 

Another assumption which we have made is that harmonics higher than the 

second degree may be neglected in the tidal disturbing function (Equation n~4). 

To show that this approximation is a good one, we note that the second degree 

. .6 

harmonics in the disturbing function are multiplied by j , where a is the 
radius of the earth and r is the earth-moon distance; this may be seen from 
Equations (II-3) and (II-4) . If third degree harmonics were included in the dis- 
turbing function, then they would be multiplied by ; likewise, fourth degree 

/a\io 

harmonics would be multiplied by (—1 , etc. Hence the third degree terms are 

reduced by a factor of frpm the second dfegree terms. “For r = 3.83a, 

/a\2 

f— j = 0.068. Also, the contribution of the third degree tides to the rate of 
change in time of the inclination is small compared to that of the second degree 
tides (see the discussion in the last paragraph of this section). Thus the restric- 
tion to the second degree terms in the disturbing function leads to only a small 
error. 

We have further assumed that the moment of inertia of the earth C had its 
present-day value of 8.11 xT0^° g-cm^ = 0.33 Ma^. This implies that the 
earth's core had already formed. Darwin assumed that the earth was homoge- 
neous (C = 0,4 Ma^), as well as incompressible, etc., for reasons of tractability 


in solving for the response of the earth to the tidal force. Changing the moment 
of inertia to its value for a homogeneous earth would lead to only slight correc- 
tions in our results. 

Heating of the earth by the dissipation associated with the friction does not 
appear to be significant. The energy deposited in the earth as the moon moves 
from the Roche limit at 2.89 earth radii to 3.83 earth radii amounts to 1.43 x 
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10^® ergs. Assuming a specific heat of 10’ ergs/g-deg, the average change in 
temperature of each gram of matter in the earth is only 24“K. 

Our most crucial assumption was that the earth behaved like a highly 
viscous liquid (v » 10 poises). Whether the earth could have behaved in this 
manner when the moon was at 3.83 earth radii depended upon the rheological 
properties of the earth at that time; these properties are imknown. 

- O'Keefe (1972) points out that since the tidal potential varies like the in- 
verse cube of distance .(Equation A-6, Appendix A), the tidal forces acting on the 
earth were 4000 times greater when the moon was at 3.83 earth radii than they 
are today, so that the material in the earth may have been near the elastic limit. 
In such circumstances the earth may have behaved like a highly viscous liquid. 

At the present time the mantle oMbe earth responds to the tidal forces like 
an. anel^tic-solid, with the tidal lag angles being small (MacDonald 1964) . How- 
ever, the mantle responds to deformations of the earth's surface caused by ice 
loads as though it had a viscosity of about 10^^ poises (Gutenberg 1959, Chap- 
ter 9), requiring thousands of years to rebound after the removal of the loads. 
(This may be ejqplainable in terms of diffusion creep; see Kauia 1968, pgs. 101- 
104). Now the period of the O tide with speed n - 2fl is given by 2 7r/(n - 2D), so 
■that the period. ranges from infinity to about 5 hours as the moon moves through 
3.83 earth radii distance. Hence if the earth has a characteristic response time 
between these two extremes, then it should be excited by the 0 tidal force as the 
moon passes through 3.83 earth radii. If the dissipation were great, then the lag 
angle of the tide would be large. Hence it is by no means clear that the earth 
would not respond as we have assumed, even with the present internal conditions 
in the earth, where the characteristic response time is thousands of years. 



60 


Our last important assumption was that the moon may have been perturbed 
out of an equatorial orbit by 2.5 to 3® at 3.83 earth radii distance from the earth, 
thus explaining the present inclination of the lunar orbit to the ecliptic. Whether 
the moon could suffer such a perturbation is not clear; conditions at that time 
may have been chaotic enough to produce it. However, several sources of the 
perturbation may be ruled out. The first obvious source of perturbation is a 
collision of a large meteoritic object with the moon. If such a collision occurred, 
then large amounts of meteoritic nickel might be expected to spatter over the 
moon's surface*. Large amounts of nickel are not observed in lunar samples. 
The third degree harmonic in the earth's figure will not give rise to long period 
perturbations in the inclination if the moon's orbit is circular- Further, it may 
be shown from Equation 38 of Kaula (1964) that the tides associated with the third 
degree harmonics in the tidal disturbing function give « 0, just as in the 
second degree harmonics, but that these terms are much less important than 
those discussed here. Also, the disturbance in the inclination caused by the pre- 
cession of the earth's axis and the moon's orbit may be shown to be quite small 
(«1°). The question of the source of the perturbation remains open. 

B. Relation of the Besults to Theories of the Moon's Origin 

We will now examine how our results relate to the theories of the origin of- 
the moon. The three principal theories, namely fission, accretion, and capture 
are reviewed by Kaula (1971) . 

Darwin (1880) proposed that a primitive body rotating with a period close to 
its natural oscillation period was disrupted into the earth and moon by resonance 
oscillations induced in it by the sun. (That this was not at all likely was shown 


*I am indebted to Dr. O 'Keefe on this point. 
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by Jeffreys 1930.) The moon would necessarily be thrown out in the equatorial 
plane of the earth. Darwin was forced to assume that the primitive earth had a 
very high viscosity to solve the inclination problem. He derived Equations {m-21) 
and (m-22) which give rate of change of inclination with distance in the limit of 
infinite viscosity. After commenting on the absurdity implied by the equations 
that the rate of change of angle was infinite when the earth rotated twice as fast 
as the moon revolved, he assumed that the viscosity merely had to be very large 
to increase the inclination of the moon's orbit to the equatorial plane from an 
infinitesimal disturbance to an appreciable angle. Darwin took this as the solu- 
tion to the inclination problem and let the matter rest. 

Our detailed analysis (Chapter III) shows that the initial perturbation in the 
inclination of the moon's orbit to the Orth's equator must be about 2.5-3.0'’ to 
explain the present inclination, with the viscosity of the earth being greater than 
10^’ poises. 

O'Keefe (1969) in his version of the fission theory suggested that the primi- 
tive body had greater mass and twice the angular momentum of the present 
earth -moon system. As the primitive body spun up, its figure progressed along 
the sequence of the well-known Jacobi ellipsoids and pear-shaped figures (Jeans 
1961) until it fissioned into the earth and moon. The system then lost mass and 
angular momentum through intense heating. While taking over Darwin's results, 
O'Keefe further suggested that even if the earth were molten after the moon and 
earth separated, the moon's orbital evolution would not begin until the earth 
cooled off appreciably, so that the moon would not arrive at 3.83 earth radii 
until the earth's viscosity was quite hi^. 

In Chapter III, Section D, we investigated the orbital evolution of the moon 
as the earth cooled off for a number of different activation energies and coeffi- 
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cients in Equation (DI-23) . In view of the wide variety of results obtained in the 
temperature of the earth when the moon arrives at 3.83 earth radii, it appears 
that the self-regulating mechanism proposed by O’Keefe does not exist. 

The accretion theory states that the moon formed from a ring of particles 
in orbit about the earth. The particles collided with each other and stuck to- 
gether, ultimately building up into the moon. 

The ring of particles would be expected to lie in the proper plane. The 
orbits of particles inclined to the proper plane would precess, thus lowering the 
chances of collision; at least all the orbits would intersect the proper plane, 
favoring accretion there. H the moon accreted from the ring much beyond 3.83 
earth radii, then the moon would tend to remain near the proper plane, so that 
its present inclination to the ecliptic could not be explained. However, if the 
moon formed injth?=p5^er plane within 3.83 earth radii (essentially in the equa- 
torial plane), then the mechanism proposed here for driving the m oon o ut of the 
earth's equatorial plane could have come into play. 

-At any rate, regardless of how the moon arrived at 3.83 earth radii in the 
equatorial plane of the earth, whether by fission or accretion, if the viscosity of 
the earth was greater than 10^’ poises, and the moon suffered a 2.5-3.0“ per- 
turbation in inclination at 3.83 earth radii, then the present inclination to the 
ecliptic could be explained. 

The capture theory has a simple answer to the inclination problem: the 
moon was captured in a highly inclined orbit to begin with and tidal friction has 
acted to decrease the inclination to its present value (Gerstenkom 1969, 
MacDonald 1964). Of course, these theories begin with the present inclination 
of the moon’s orbit to the ecliptic and solve the equations of tidal friction back-' 
ward in time to discover the moon’s inclination at capture. 
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These theories once again assume weak tidal friction with small lag angles. 
However, if the viscosity of the earth is greater than 10^^ poises, and the moon 
arrives at a distance less than 3.83 earth radii in an inclined orbit, then the 
inclination must drop below the critical angle before the orbit can expand 
past 8.83 earth radii (Chapter HE, Section E). Thus for large viscosities 
(greater than 10^^ poises) the orbit becomes nearly equatorial and we are faced 
with the same problem as before. 

C. Summary of the Important Results 

Assuming that the earth behaves like a viscous liquid in responding to the 
tidal force, and that the moon is in a circular orbit about the earth, with the in- 
clination of the lunar orbit to the earth's equator ~20®, then:* 

(a) If the moon is less than 3.83 earth radii distance from the earth and 
the inclination of the orbit to the earth's equator is steep, then the orbit 
may contract and then expand, provided the viscosity of the earth is 
greater than 10^’ poises. The orbit will expand monotonically if the 
viscosity is less than 10^® poises regardless of the inclination. 

(b) The inclination of the lunar orbit to the earth's equator will decrease 
or remain zero if the moon is closer to the earth than 3.83 earth radii, 
regardless of the viscosity. 

(c) The inclination of the lunar orbit to the earth's equator must be less 
than 2.7° when the moon is at 3.83 earth radii if the earth's viscosity is 
10^® poises; at higher viscosities the inclination must be even lower. 

(d) The lunar orbit will expand monotonically if the moon is at a distance 
greater than 3.83 earth radii from the earth, regardless of the viscosity. 

*<20* so that Equations (Ul-9| through (lll*t2) hold. 



64 


(e) The inclination of the lunar orbit to the earth's equator will increase, 
or decrease (or remain zero) for earth-moon distances greater than 
3.83 earth radii depending upon whether the viscosity of the earth is 
greater than, or less than 10^^ poises. 

(f) If the viscosity of the earth is greater than or equal to 10^® poises, and 
the plane of the lunar orbit is perturbed about 2.5 to 3 degrees out of 
the equatorial plane of the earth when the moon is just beyond 3.83 earth 
radii, then the present five degree inclination of the moon to the ecliptic 
may be explained. 



III. Comments on Atmospheric Regractton 


In order to generally address the questions on the formation of the 
earth-moon system which were raised in the analysis presented by the previous 
section, precise measurements of the current motion of the earth-moon system 
are required. The Principal Investigatory has been involved with ^several 
of the current techniques which are either now in use or are candidates 
for future application (i.e. Laser Ranging to Lunar Corner Reflectors 
Laser Ranging to Artificial Satellites, classical Astrometric techniques 
with PZT^s and Very Long Baseline Amplitude Interferometry in the optical 
domain) . 

All of these techniques are, or will be, limited in their ultimate 
.accuracy by large scale refractive phenomena in the earth’s atmosphere. 

Prior to this grant, a preliminary study was initiated on special techniques 
to use the optical dispersion of the earth^s atmosphere, in conjunction 
with measurements performed at more than one wavelength, to sense the 
magnitude of this large scale errors and provide data for their correction. 

Under this grant, these general methods were studied in more detail. 
Certain general problems on calibration of the water vapor content were 
also addressed. 

As an illustration of this type of work two subsequent studies, to 
^which the studied conducted under NCR 21-002-390 contributed and which 
used the tools developed under this grant, are attached as appendices. 

The first appendix contains a theoretical analysis of the techniques 
of using observations at Two Color Technique for measuring of the atmospheric 
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refract ion ♦ This data has the format which may be used to correct a PZT 
observation for this large scale effort* Appendix I also contains a dis- 
cussion of a hardware technique to implement this theoretical analysis* 

The second appendix discusses the use of the basic Two Color 
Techniques (by which one may use the measurements at two different wave- 
lengths) in Very Long Baseline Amplitude Interferometry* Here again this 
data may be used to correct the systematic error introduced by the 
atmosphere. This section generally reviews the technique. This portion 
of the work under this grant (NGR 21-002-390) , consisted of the exploration 
of a number of the important specific questions which are raised in this 
discussion. 

In addition, several analyses were conducted concerning empirical 
techniques to correct for the systematic atmospheric effects which are 
related to the lunar motion, as they applied to laser ranging application. 
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I. INTROnuCTTON 

Atmospheric refraction^ with Its non-prcnUrniblo vnrffUfnnH» is ono 
of the dominant sources of error in the precise de(erminntion of the 
apparent position of a star* It appears that this will be t he 'dominant 
error source for the new generation of instrumentation now being Installed 
for the detorminatiou of stellar catalog and Universal Time (l.e*^ tlie 
USNO 65**em PZT). The e.ffeets of anotnolous refraction measurements of 
stellar position propogate to the various quantities derived from this 
data. This includes tlie compilation of stellar catalogs, the determination 
of Universal Time, and variations in the latitude, 

A unique instrument, the Two-Color Rcfractometer, can, for the first 
time, directly measure the total deviation which the light of a star 
suffers due to refraction and thus can be used in a measurement program to 
evaluate the various aspects of anomolous refraction- In addition, this 
system may be used to make fundamental astrometric observations which are 
free from the deleterious effects of refraction- 
A- Atmospheric Refraction 

There are two aspects of refraction errors which are useful to discuss 
separately- The first of these effects are the short-term errors or the 
"random error** in an individual measurement. These errors affect an individ 
ual measurement so that it may be significantly different from measures 
made before and after it on the same night- This type of error may be of 
the order of 0-15 arc-seconds or larger- When we later discuss all of the 
errors more generally in terms of the power spectra of the anomolous zenith 
refraction, the random error is characterized by the magnitude of the power 
spectra in the domain having periods of a few minutes. 

The many individual measurements may be analyzed and the results 
averaged together for an entire night. This procedure should result in 
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an improved determination of the quantities derived from the stellar 
position measurements • The improvement should be parameterized by a 
factor proportional to l//n* This will occur if the errors of each 
separate measurement are independent. However, if there are phenomena 
which affect the anomolous refraction which have long periods, then the 
determination of the averaged quantities may now show this improvement. 

Again, discussing the question in terms of the power spectra, one will 
obtain the improvement in performance only if the power spectra vanishes 
for periods longer than a few minutes. We now address the effects of the 
power spectra for periods greater than the time which is required for a 
single observation. The systematic error to be expected may be derived 
from the power spectra and it will be dominated by the magnitude of the 
power spectra with the period in question. 

B. Derivation of Star Catalogs 

We wish now to consider the question of the improvement of the star ■ 
catalogs. This will require one to both improve the measurement of the 
single observation and reduce the systematic errors obtained when averaging 
data over T-ong periods. The former will improve the basic random error 
and the latter will insure that the result of many measurements will reduce 
this diminished random error in proportion to l//n. In general, the stellar 
catalog measurement using the Transit Circle is most affected by three types, 
of error. These are: the anomolous refraction, the encoder errors, and 

the lack of stability of the personal equations. It is difficult to 
separate the effects of these errors prior to the measurements of anomolous 
refraction. Therefore, we shall address the remaining remarks to measure- 
ments made with the Photographic Zenith Tube (PZT) . Por the PZT, the 
errors of encoders and lack of reproducibility of personal equation are far 
less significant than in the Transit Circle. In addition, to first 
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order one is independent of normal refraction and only sensitive to the 
anomolous refraction* In this discussion, we refer to the “normal refrac- 
tion“ as that theoretical quantity derived from the zenith distance of 
the star, and the local pressure, temperature, and humidity* The anomolous 
refraction is the difference between this normal or calculated refraction 
and the actual measured refraction. 

In the case of the PZT, the normal random error is of the order of 
0.18 arc^seconds. This is the internal precision for one night and 
does not include the effects of longer period phenomena which is probably 
significantly smaller. It is more difficult to estimate the errors with 
a longer period. These effects, due to phenomena like the heat island 
effect of Washington, D.C., the seasonal variation of the heat island 
effect, and long term weather effects may lead to violation of the 
horizontal uniformity which is assumed in the normal data reduction 
procedure. The derivation of the averaged measures for one night at 
the USNO may be of the order of 0.06 arc'-seconds for periods between one 
and one hundred days. 

C. Two-Color Refractometer 

In order to address the role of anomolous refraction in these 
determinations of star position, we suggest the development and use of 
the Two-Color Refractometer (TGR) . The TCR is a system which measures 
the magnitude of the refraction due to the atmosphere when observing 
a star. This measurement is not affected by problems of an encoder, 
a star catalog, or variations in personal equation. It is a measurement 


which isolates the refraction alone. 
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The primary subsystem of the Two-Color Refractometer is the Quadrant 
Sensor System (QSS) . The QSS is an operating system which has been 
developed at the University of Maryland* The QSS is an “eyepiece” which 
permits one to determine the apparent star position with great accuracy. 

This determination of the apparent direction to the star is performed by 
observing in two different colors. Using a knowledge of the index of the 
refraction and the dispersion of the air, one may compute from this data 
the instantaneous magnitude of the angle of refraction. Using the local 
temperature, pressure, and humidity, one may then compute the normal 
refraction. The difference between the normal angle of refraction and the 
measured angle of refraction is the anomolous refraction# The latter will 
be done by exactly the same calculation which is normally used by USNO 
for the reduction of transit circle, i)bservations. 

In general, the anomolous refraction will depend upon both zenith 
distance and time. Let us for the moment neglect the time variation. 

Then a number of stars may be observed "at one instant” (actually a 
period of about thirty minutes) . The measured value of the anomolous 
refraction may then be expressed in a power series function about the zenith 
in terms of "east” and "north" direction. The first term, which will be 
a constant term, which might be due to the use of erroneous values of 
local pressure, temperature and humidity, for the computation of the normal 
refractions. Any linear terms which appear would be due to gradients in 
the atmosphere. Quadratic terms should be small if the theory is proper 
and the constant term is not large. 
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II. OPERATION AND PEEPORMANCE OF A QUADRANT SENSOR SYSTEM 
In this section we consider an ideal Quadrant Sensor System (QSS) 
and discuss the expected performance of such a system. The objective of 
the first part of this section is to define the type of sensor system 
which would be most desirable and to provide a model for the performance 
analysis which will appear later in this .section. The realization of 
the QSS system will be discussed in Section IV. 

The QSS is based upon a special photosensor, the "Quadrant Photo 
Sensor" produced by the Electronic Vision Company, a division of Science 
Applications, Inc. The Quadrant Photo Sensor (QPS) consists of a quadrant 
silicon diode which is mounted in an evacuated envelope. Photons are 
converted to photoelectrons at the photocathode and accelerated electro- 
statically to an energy of 15 KeV. They are then electrostatically imaged 
onto the quadrant diode. The overall mechanical configuration of the 
Quadrant Photo Sensor (QPS) is illustrated in Figure 1. 



Schematic Diagram of Quadrant Photo Sensor 
Figure 1 

The photoelectron, which is accelerated to 15 KV, produces extensive 
ionization within the silicon in the photodiode. The output of a photo- 
diode which is bombarded by a single photoelectron is a charge packet which 
consists, at the normal operating voltage, of about 2,000 electrons. This 
pulse signal generated by this charge packet is amplified, shaped, and 
finally detected with a level discriminator. In order to operate such a 
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system, we wish to provide functions indicated in Figure 2. 



QSS Detection Subsystem for' Single Photoelectrons 
Two of the Four Channels are Shown 
Figure 2 

Thus Figure 2 illustrates the photon detection system, which provides a 
high-level pulse for the arrival of each photoelectron. By accumulating 
these counts, one would have available the intergal number of counts for 
any interval. However, we wish to use such a system in several applications 
where a greater variety of outputs are required. The above system will 
be used as the basis for the calculations presented in Section II. At 
present, however, we shall now continue this discussion to describe the 
full "Quadrant Sensor System requirements. The subsystem described in 
Figure 2 is connected to a data processing subsystem. 

The QSS has been designed both to provide the input for a telescope 
used for the tracking of stars, but also to detect laser returns from 
satellites and track on these returns. However, the latter function is 
not used in this application and will not further be considered. 

When operating in the photon sensing mode, the system is expected to 
produce an output signal which is proportional either to the Intergal number 
of counts or the count rate smoothed in some fashion. Two outputs of this 
type, which are tihe difference between the counts from each axis, are 
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available (X high pass, Y high pass). In addition, the output is then 

filtered with a low pass filter having a time constant of about 1 second. 

This is available as X , and Y , , The high speed outputs 

low pass low pass 

will be used as the input for an external subsystem to -stabilize the 
image with a high frequency response. The low frequency output will 
provide the signals required to move the telescope. The signals described 
above are available in analog form in order to operate the analog servo- 
loops. They are also multiplexed onto an A-to-D Converter to provide 
outputs which may be used for computer control and recording. Thus this 
unit is capable of both stabilizing the image, and through the photon 
counts, provide interval number of counts in the various channels. 

In this discussion we shall evaluate some of the numerical performance 
parameters for the quadrant sensor. In this section we shall presume 
the use of a quadrant photon counting sensor. We will assume that the 
photocathode has a uniform quantum efficiency and that the electronics 
system imposes no rate limits. The description of the actual sensor and 
its limitations appears in a later section. 

Photon Noise in Direct Mode 

In this section we shall compute the photon noise which affects the 
determination of the centroid of an image by a Quadrant Sensing System. 

For pedigogical ease, we shall first consider a stellar image which is 
square, so the brightness has the form indicated in Figure 3. 



Brightness of Stellar Image 
Figure 3 
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The difference between the performance criteria for this approximation 
to the stellar image (by a square pillbox) and a more realistic approxima- 
tion, l.e., a Gaussian or a Lorentzian form, will be discussed later* 

The more realistic forms will result in about a 2% modification of the 
numerical results obtained for the square pillbox image. 

The results shall be expressed in terms of the total number of 
photoelectrons which arrive on all four quadrants (denoted by the symbol 
N) during a time Interval denoted by the S3rmbol T, The full width of the 
image described in Figure 3 is denoted d, which might typically have a value 
of two arc-seconds for a reasonable telescope used at a good astronomical 
site. The one dimensional projection of such an image has the form 
indicated in Figure 4. 



One Dimensional Projection of Square Pillbox Image 

Figure 4 

The number of photoelectrons per linear arc-second at the center of the 

. f\j 

square xmage, denoted by N, is given by 

Sf = N/d = N/(d /3) 
o 

where is the ”rms diameter*’. This quantity is obtained by doubling the 
rms radius, which is the deviation of the light Intensity distribution over 
the image plane. We may reexpress N in the form 

N = (n N)/(2d ) 

0 
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where we. have defined the symbol n to represent a numerical coefficient 
which is related to the shape of the image. Thus from the above expressions, 
we see that for a square image 

p = 2//3 = 1.547 
s 

Now we consider an image with a more realistic profile, i.e., a 
Guassian. The one dimension projection of such an image is indicated in 
Figure 5. 







0n£ T)mmiiyncii'Prtfj£ctm ^ Gaussian Jmaje 


Figure 5 

The intensity li, expressed as the number of photoelectrons/arc-seconds, has 
the form 


Thus at origin 


NCx) 


e 

d 

o 


<v <v 
N = N(o) 


/2 N 
vV d 


= n N/2d 
g o 


Where we have again used the coefficient as p. Thus for the Guassian 
we have 
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Thus in general, the value of n will be dictated by the tnodel for 

the seeing disk, or from a more pragmatic veiwpoint, from an actual 

measurement of the size and profile of the seeing disk* 

We may later take into account a large deflection or error dn image 

position by expressing n as a function of the offset of the image. Thus 

if there is an average offset of A , then n (A) will have a smaller value 

Xn this x^ay we will be able to address some problems which will arise in 
the two color application. Thus for a Guassian image 


„ /.V _ 2/2 ~ 2 A^/d ^ 

n (A) e 0 




The quantity N is physically related to the change in signal or differen- 
tial count rate for motion of the image. Thus if we apply a sinusoidal off-- 
set or a "dither^* to hhe position by an optical device, we will see a 

r\j 

similar variation in N* This will obviously be related to the profile, 
as expressed by n# On the other hand, if we assume a profile, we may 
invert this relation and use the measured value of N which is obtained 
by applying a dither to the servo-system to evaluate the seeing disk 
diameter, using the expression 


d = n N/2N 
o 

If the image is nominally centered, then a small image offset which 
is parameterized by the value of 0 will result in a difference in count 
rates which is given by; 

M = Nl~ + H-0) - (Nj,^ - N-0) 


We will describe the one-dimensional case. If the image was initially 
centered and the quantum efficiencies of the effective areas of the photo- 
cathode are equal then N = N and the above equation becomes 

IjO ko 



which may be reexpressed as 


AN/N - n 0/d 
o 

or using AN to compute the offset angle 

0 » (AN/N) d^/p 

We may now define a relative angle which is expressed in terms of the image 
diameter, so we have 

a = 0/d = AN/Np 
o 

This expresses the angular offset due to an inequality in counts as a 
dimensionless angle - 

Now let us determine the standard deviation of the angular position 
due to shot noise, i,e., the statistical error in the determination of the 
centroid due to the photon statistics. In each pair of quadrants (presuming 
we sum the photo counts from both quadrants of each side to determine the 
error) the mean value is 

Thus for large values of N, the statistical noise in the value of each of 
these is 

6Nl = = Ml 

Since the statistical noise in each of these quantities are independent, 
we have 


5N = 
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inserting this value of 5N into the above expression for a yields an 
equivalent angular offset due to the noise. Thus we have 

a ■= l/n 

At this point, these considerations will be extended to include the 
effect of a uniform sky background on the performance. Thus the presumed 
brightness now takes the form 



'Distn'^iion jor 

VKik ‘Vkctosensdwe N\arkei 

Figure 6 

where N, is the total nimber of counts which appear in the four quadrants 

D 

due to the brightness of the background sky. This is equal to the product 

\ 

of the field of view (in arc-seconds) , the background illumination (in 
photoelectrons/arc-second/sec), and the integration time (in seconds). 

The product is expressed in photoelectrons. The image of the star is 
characterized as in the previous discussion. In this treatment, one 
simply adds the background to the image,, so the mean number of counts on 
each side is given by 

Nj, - Nj = (N + \)n 

Thus the photon shot noise in one of the pairs of quadrants (again for 
relatively large signals) is given by: 
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6N. = == /(N + Nj/2 

L L D 

and the values for and have similar expressions. Now if we take the dif- 
ference of these two Independent qualities and use this quantity to determine the 
apparent error in pointing angles, with the presumption that the photon 
statistics on the two sides are independent and the sky Illumination is 
spacially uniform, we have 


AN = Al + N, 
0 


Thus for a we have 


a = /N + N, 


VnN = 


/i~+ 


(N^/N) 
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where we have used the total counts for the noise in the numeration, and 
the counts in the image for the signal which appears in the denominator/ 
We have two types of behavior of interest, the first of which consists 
of the case when the background is negligible. For this case the error 
becomes 


a - l/n 


as discussed earlier. This is illustrated for the case of a Gaussian image 
in Figure 5 by the solid line. 
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The other extreme case is when there is a l-high level of background 
illumination. This might represent observations conducted in the daytime 
with a blue sky. In this case we have: 

a = /N^/N ^ j nN 

This may be summarized in graphical form by the dashed line in Pigure 7. 

Thus Figure 7 Indicates the performance that one may expect to obtain 
with the quadrant sensor, fhe expected error in the relative angle and 
the expected error in arc-seconds, presuming a seeing disc diameter of two 
arc-seconds are related to the total nutdber"of photo electrons received 
during the observation period. 

We may now transform this relationship to describe the behavior with 
respect to stellar brightness and telescope aperture. Let us, for example, 
consider a 48-inch telescope and photoelectron count rates which are caused 
by a particular stellar brightness for a G-2 star (or solar reflection 
from a satellite) . We assume a quasi- S-20 photocathode as measured for 
a tube of this type and six mirrors with reflectivities of 0.7, 0.7, 0.86, 
0.86, 0.86, and 0.86. Thus for the case of a Guassian image with negligible 
sky background we have the relation indicated in Figure 8. 
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This indicates the ability to go to rather faint objects and still 
have very small angular errors due to the photon counting. The numbers 
in parentheses indicate the results for a 4.8-inch aperture. 
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III. TWO COLOR APPLICATIONS 

.The Quadrant Sensor System, as a component of the Two-Color Refrac- 
tometer system may be used to obtain high accuracy measurements in 
absolute astrometry. In particular, it will determine the deviation in 
the apparent star position caused by the earth's atmosphere. As star 
light enters the earth's atmosphere, the deviation of index of refraction 
of the atmosphere from that of the vacuum results in a change in the direc- 
tion of propogation. For a flat, horizontally-stratified atmosphere, 
this change in direction is a known function of the index of refraction of 
the air in the immediate vicinity of the telescope (the local index of 
refraction) . This may be determined from the local temperature, pressure 
and humidity. However, if there are horizontal inhomogeneities, this 
change of propogation direction can not be predicted from a measurement of 
local surface conditions. In this section, we address a procedure to 
determine the magnitude of this change in direction (the "angle of refraction") 
from measurements made on the light of the star using the Quadrant Sensor 
.System. 

We take advantage of the fact that the magnitude of the angle of 
refraction depends upon wavelength, due to the dispersion of the index of 
refraction of air. Thus hy measuring the apparent position of the star 
as viewed in red and in blue light, the difference in apparent position 
may be evaluated. This dispersion angle may be used with a knowledge of 
the dispersion in the index of refraction to determine the angle of refrac- 
tion. Thus using the dispersion angle between the red and blue images, 
and the known dispersion of the air, we may determine the 'deviation of the 
direction of propogation of the starlight due to the atmosphere. 
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In order to illustrate the operation of the Two-Color Refractometer 

(TCR) let us consider the basic relations governing the refraction. The 

symbol 0^^ describes the deviation of a ray of wavelength X from its 

initial direction and (f> is the true zenith distance of the star. The 

index of refraction of the air at the wavelength X is denoted by n . 

X 

Thus from Snell's law, we have 

= ("x 

, A 

Since is very close to unity, -we shall approximate this expression by: 
where the definition of these quantities may be seen in the following figure 



B^ecf c^'K^acfwn c*i 

9 . 


While the figure illustrates a simple idealized case, a little 
reflection indicates that equations apply to the more general case of a 
continuously variable, non stratified, spherical atmosphere. For the 
calculation made in the later part of this section, we use the following 
expression for the index of refraction for dry air 

•" -2 1 

(n^ -1) X 10"^ = 64.328 + 29498.1 (146 - Xq ) + 255.4 (41 - Xq ) 
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This applies for a temperature of 15”C, a pressure of 760 mm and- no water 
vapor. This expression is not the optimal form for expressing the index 
of refraction when reducing the data for the operational system, but it 
is entirely satisfactory for the system analysis and error analysis 
presented in this section. 

In order to permit the independent variation of the wavelength selected 
for observation within a reasonable error analysis, we shall use the blue 
and red measurements and a parameter which depends on the wavelengths used to 
determine the magnitude for the angle of refraction at a reference wavelength 

O 

chosen to be the center of the visual band (5600A) . Thus for the wavelengths 
at which we will measure the relative angular position, we use: 

0R = -1)^> 

We will then use these expressions to determine the zenith distance ip, 
which we will relate to the angle of refraction in the green by the relation 

Thus we will determine the angle of refraction in the visual band. This 
particular equation does not appear to represent the optimal procedure 
for minimizing thd propogation of experimental errors but we shall use it 
as a useful tool for system analysis. Thus we obtain an expression of the 
form: 


9, -MSb-Sr) 

where 


= (n^-D/ECng -i)- - (dr -D1 
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If we measure the difference in apparent angular position of the blue 
and the red image, this angular difference, multiplied by the quantity Cl, 
will yield the angle of refraction suffered by the visible light. This 
is the value for local conditions which .are 15“C and a pressure of 760 mm Hg. 

In order to provide a general frame of reference, let us consider 
the numerical values for the refraction and the dispersion angles. 

Let us consider the value of K for various choices of the limiting 
wavelengths. We have 


Effective Blue 

Effective Red 


Wavelength 

Wavelength 

Q 

3500A 

6000A 

30.36 

• 

3300A 

6000A . 

25.27 

• 

3500A 

7000A 

26.89 


Relation Between Selected Wavelengths and n Multiplier 

Table -1 

Since the eventual error is proportional to the value of n, we would 
prefer the smallest possible value. From Table 1 we see that the value 
of Q decreases and thus the error multiplier decreases as the spread 
between the wavelengths Increase. This dependence is particularly 

O 

sensitive to the wavelength of the blue observation where a change of 200A.-. 

O 

pctrvides as much improvement in Q as a change of lOOOA in the wavelength 
of the red observation. ' However, for various technical reasons we will 
aceapt the first case as a reasonable, practical limit. 

Let us now evaluate the sources of error in this determination. 
Differentiating the expression for 0^, -we have: 

®V-f 
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We now presume the uncertainty in the measurement of the angular 
position of .the image in the red and in the blue have the same magnitude, 
and are independent. Concerning the former, it assumes similar filter 
widths which will probably not be the case, but It is a useful point for 
discussion. Thus for the part of the variation due to uncertainty in 
angular measurement we have 

50y vT n 60 



We may now use the relations between the total number of counts and angular 
errors to express the accuracy of the determination of 60^ as a function 
of the total number of photoelectrons. In Figure 10 we have also indicated 
the various system operating points obtained when the QSS is running at the 
maximum count rate (the normal procedure) for various intervals of time. 
These results include the effect of the UV filter and the red filter. 



iciat iiumher (f fheietis 

'RefaUnt hituttn ^oUl timler ^ cpunh and error tngnmoim nfrorim fr TCR. 


Figure 10 
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We now consider the other source of error which appeared in the 
Equation for 60^. This is the possible variation in the quantity It 
is presumed that ^ is calculated for the local value of refraction (i.e. 
using pressure and temperature) . The predominant remaining reason for 
an uncontrolled uncertainty in 60^ is due to the variable water vapor 
content in the atmosphere. In order to evaluate this phenomenon, the 
following table is presented. This presents the value of 0, and the change 
produced in Si due to precipitable *water for various pairs of effective 
wavelengths. We use the mean water vapor which is expressed in the terms 
of precipitable water vapor (for a flat atmosphere) . 


Effective 
Wavelength 
for Blue 
Filter 

Effective 
Wavelength 
for Red 
Filter 

^ for no 
Water 

Precip- 

itabie 

H^O 

SI with 

Water 

Vapor 


3500 

6000 

30.3657 

10 mm 

30.1885 

0.0058 

3300 

6000 

25.2696 

10 mm 

! 

• 25.1232 

0.0058 

3500 

7000 

26.8927 

10 mm 

26.7351 

0.0059 

3500 

6000 

'30.3479 

1 mm 


0.00065 


The Effect of Water Vapor on Systematic Errors 

Table 2 

Thus we see that the water vapor correction is approximately proportional to 
the amount of water and generally independent of the selection of the effec*- 
tive wavelength of the filters. 

Let us now consider a numerical example. If we observe a star at 
zenith distance of 30* as one might with an astrolabe, then using the 
expression for the refraction under standard temperature and pressure 


= 58V2 tan Z 
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which applies to a flat atmosphere we have, at30t a deflection of about 34". 
If we wish to be able to determine this correction with an accuracy of 0V05, 
then we need 


69 

V 


0V 


60 , 0V05 

SS I ■ ■ I I 

ft 33V6 


0.0015 


This means one to know the water content to 2.3 mm of precipitable water. 

However, this requirement is not particularly difficult since 
commercially available microwave radiometers arc accurate to about 
1 or 2 mm of precipitable water. In addition, a system ttached to the 
TCR using near IR measurements in and out of a water absorb t ion band 
should be significantly more accurate and less complicated. 

Some Practical Considerations 

We now consider a few practical aspects of the design of the TCR. The 
full discussion of the system will be considered at a later time. 

We first consider the procedure for sampling the light at two different 
colors. One might attempt to measure the position of the red and blue images 
at the same time with two different quadrant sensors. The parameters which 
enter Figure 6 are based upon this assumption. However, a more practical 
approach to this problem is to use a single quadrant sensor. The detailed 
reasons will be considered in a separate discussion. We shall, briefly, 
discuss this approach and how we should use it to maximize the accuracy of 
the evaluation. Thus we note that due to the fact that the refraction in the 
atmosphere has significant power at various frequencies, one might wish to sample 
the two different colors continuously and at the same time. If this is done, 
then both the power at high frequency and the power at low frequency (the image 
motion) are properly tracked and do not affect the accuracy of the measurement. 
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However, the mechanical stability requirements strongly suggest the 
use of a single detector. If one switched the two colors onto the one sensor 
at a rate of several hundred times a second, one may show that the residual 
error is expected to be random and has ,a value of less than 0.01 arc- 
seconds for a measurement lasting several minutes. 

The remaining uncertainty is in the determination of the water vapor. 
From the previous discussion, and typical values of water vapor to be 
expected, this effect may be totally ignored for the initial measurements 
out to a distance of 15® from the zenith. 

Expected Accuracy for Measurements of Zenith Refraction 

We now discuss the accuracy which we would expect to obtain with 
measurements using this instrument. Thfese xd.ll be expressed in terms of the 
measurement errors to be expected on the power spectral density, and xd.ll 
be based upon use of the 48-inch telescope at the Goddard Optical Research 
Facility (GORF) . However, approximately the same numbers apply to the set 
of measurements which would be performed on the 36-inch telescope. For 
these evaluations, we will presume the use of the present circuit for 
the quadrant which, has a raaximxjm counting rate of 10 kilohertz/ 

quadrant. This will permit the observation of 8th magnitude stars with 
the proper count rate. Since we will limit our measurements to a region 
that is between 5 and 10 degrees, it is worth noting that we would 
statistically expect to have about 80 stars within a five degree circle 
of the zenith. 

In general, the relationship between the stellar magnitude, the tele- 
scope aperture and the expected photon statistical error for a ten minute 
observation is expressed in Figure 11. 
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Figure 11 

Combining the results of Figure 10 and Figure 11, we see that we may trade 
off between the aperture size of the telescope with the observing time to 
use a small telescope and still have a reasonable number of stars to observe. 

Now presuming this count rate, we can evaluate the power spectrum of 
the system errors which we expect. To facilitate this present discussion, 
we will not consider the true power spectrum but indicate the standard 
deviation of the system error for a measurement of a given length. The 
overall results are thus summarized in Figure 12. The system error caused 
by the photon statistical noise is illustrated by the broad solid line. 

For periods shorter than one second, one may increase the centroid of the 
white light image and rely on the smoothness of the telescope tracking 
to determine the refraction (called image in this domain) rather than the 
two-color extraction. This is expected, with some study of the records 
and knowledge of the telescope, to give about 0V2. If these measurements 
were conducted on the 100-inch telescope at Mt, Wilson, which has a very 
smooth and regular drive, and were conducted on a night of low wind, this 
system error could probably be reduced to 0.1 arc-seconds. For periods of 
more than 10 seconds, however, one would expect that errors in the telescope 
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drive will become too large. In order to study the behavior for longer 
characteristic times, we consider a limit for a single observation of twenty 
minutes which results in an error of 0.01 arc-second. Thus one might expect 
to obtain an increase in accuracy by about a factor of three, but this will 
mean the observation of new stars. The feasible accuracy in this region 
must be more carefully explored. For the longer period term, one averages 
the data obtained on successive nights to obtain these results. However, 
since we do not expect to observe for 24 hours a day, there is a break in 
this curve as one goes to the dashed portion. Prestiming an observation 
period of two hours for a night, we obtain the dashed curve. At about 0.03 
arc-seconds one expects to run into systematic errors. It seems possible 
that these may be calibrated to the 0.001 or 0.002 arc-second level. 
Accuracies beyond this cannot currently be projected for the first genera- 
tion instrument. 



Jni^raHfin (secot^i) 


Figure 12 
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TV. IMPLEMENTATION OF THE QUADRANT SENSOR SYSTEM 

In this section we describe the Quadrant Sensor System which has been 
developed within the Amplitude Interferometry Program at the University 
of Maryland. This system is the practical realization of the hypothetical 
QSS which was described in Section TI. The basic element of the system 
consists of two standardized units, the photo-sensor head and the control 
unit which have been used in three different astronomical systems. The 
quadrant sensor system is basically a special system which detects light 
in four channels. These are closely placed and are mechanically and 
geometrically rigid. In order to satisfy the system requirements, a 
special device has been obtained from the Electronic Vision Company and an 
electronic system has been developed at the University of Maryland to make 
these observations. The requirements on this optical/electronic system 
demand a very high precision, electronically, optically, and mechanically, 
as well as a considerable amount of generality. The electronic system 
will permit the quadrant sensor to interface "to a variety of data processing 
systems, i. .e, computers or hardwired processors. 

The Photosensor head contains the Quadrant Photisil. The structure 
has already been Indicated in Figure 4 and is shown here 



Connection to Quadrant 
Diode 

Quadrant Silicon Diode 


Ceramic Header 


Schematic Diagram of Quadrant Photo Sensor 
Figure 13 
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As discussed in the first section, the circuitry following the Quadrant 
Phot os il is indicated in Figure 14. 



QSS Detection System for Single Photoelectrons 
Two of the Four Channels are Shown 
Figure 14 

\ 

This results in a pulse height distribution which permits a straight- 
forward discriminator. The preamplifier is similar to that photon counting 
Digicon system developed by E. Beaver at UCSD. The pulse height distribu- 
tion at 16 KV is shown in the following figure. 
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Figure 15 
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This system is expected to operate In environments with high levels 
of external EMI. This is a serious problem since the signals are smaller 
than conventional PMT's by a factor of about 1000 to 10,000. 

Special shielding and grounding procedures are required in order to 
prevent external TVI from interfering with the discrimination process. 

This has been provided for the input high voltage leads by the use of 
solid aluminum shields. The amplifiers are specially isolated and 
grounded. Even so, the pickup is noticable under certain circumstances. 

The preamplifier in the QSS head provides a signal which is rela- 
tively immume to external problems. This is then sent by four short cables 
to the QSS Control Unit. Here the signal is reamplified and sent to 
discriminators. These discriminators may be adjusted independently for 
each of the four channels. These are normally set at about the point 
indicated by the arrow in Figure 15. The shape of this pulse height 
distribution leads to a relative inmunity from effects due to the varia- 
tion one expects in operating voltages and circuits in the field. 

In order to interface to a computer and provide data for the closure 

of a servo-loop, the four signals are multiplexed and then applied to an 
A-to-D converter. The output is an eight-bit word. The operation acts 
as an asynchronous peripheral which supplies data when asked. Thus the out- 
put of the ADC is available as a digital signal. This generalized system 
has been used in several applications. We Shall now describe how it has 
already been used. 

Application in Single Aperture Amplitude Interferometer 

The first application of the quadrant sensor system was to provide 
Image stabilization in Single Aperture Amplitude Interferometry (SAAI) . 

In this system it operates in only one axis. It incorporates a device 
which provides high frequency stabilization with a frequency response up 
to about 60 hertz. The gain is set to permit the utilization of the 
stabilization to this level. The expected noise performance was discussed 
in the previous section* 
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As an indication of the proper performance of this system in Amplitude 
Interferometry, the measurements of the fringe visibility which are per-' 
formed when this system is operating are significantly improved* Special 
procedures have been developed, both for the proper alignment of this 
system and to permit the simultaneous use of the manual guiding of the 
telescope* This is required at present since the QSS has not yet been 
interfaced to the large telescopes due to the necessity of performing the 
coordinate transformation without a computer* 

48"lnch Telescope at the Goddard Space Flight Center 

The Quadrant Sensor System (QSS) has been used to- make an Automatic 
Guider System (AGS) to control the 48- inch telescope at the GOBI* In 
'this case, the output from the high frequency error signal is multiplexed 
and digitized and fed directly to the computer. There this data on the 
error in pointing is processed and filtered in a digital form* This 
system has succeeded in providing precise pointing information and also 
has been useful in detecting difficulties in the encoder system. Insta- 
bilities in the computer processing of the data and the telescope setup 
drive can cause oscillation under some circumstances* However, the 
computer procedure for handling the data is now being reprogrammed to 
handle this properly. 

In addition, the GSFC system also provides for the capability to sense 
the offset of lasers from a satellite reflection. To this end, the 
capability is provided to gate the photo cathode off while the outgoing laser 
pulse is near the telescope with a high voltage pulse applied to a grid 
of a modified Quadrant Photosensor* The sensor system can detect the 
arrival of a laser return* The signal on each diode is amplified, and 
then compressed for greater dynamic range, and stores before it is presented 
to the multiplexer. This sequence is activated by the receipt of a photo 



pulse which has a value which is greater than three photoelectrons. 

Thus this gives the error dn pointing with respect to the laser return. 
This system has been implemented and tested in a laboratory but the 
required computer programs for the 48-inch telescope and the required 
telescope time are not yet available in order to test the laser tracking 
system on the telescope. 
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ON THE USE OF WIDEBAND OPTICAL VERY LONG BASELINE INTERFEROMETER 
FOR THE DETERMINATION OF THE POSITION OF THE POLE 



On the Use of 


Wideband Optical Very Long Base Line Interferometer 

for the 

Determination of the Position of the Pole 
Local Reference Frame 

The discussion of the stellar measurement and the atmospheric 
correction have related 'the reference points on the piers to the 
stellar frame. .Let us now consider the relation of these two points 
to the earth. Thej^ are set into concrete piers which are buried in 
the mountain. The questions concerning these reference points are; 

L) The piers may move with respect to the local 
g^lo-giSal formations. 

2) The local formation may move (distortion of the.shape 
of the mountain). 

3) The mountain may move with respect to surrounding terrain. 

4) Continental drift may be moving the observatory, 

5) Polar motion changes the declination. 

These possible motions are all "problems" but they also repre- 
sent phenomena of geophysical interest.' The primary difficulty is 
when all these motions occur simultaneously, so that it is difficult 
to evaluate them separately from the declination measurements. To 
this end, three methods may be used to relate the reference points 
to the earth. These are: 

A) Piers — the piers tie the reference points to the local 
geological structure. 

B) Local Gravity Vector — the reference points (and the 
declination angle) may be tied to the local gravity 
vector. 
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C) Distant References — the reference points may be related 
to piers over a distance of about 10 kilometers, forming 
a net o.f local references. 

Thus the large scale motion may be separated from possible local 
motion and changes in tlie direction of the loccil gravity vector. 


Douglas G, Currie 
30 December 1971 . 
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Betermination of the Position of the Pole 
Atmospheric Correction by the Use of Dispersion 

The following is a description of a procedure to eliminate the 
effects of the atmosphere as a source of systematic error in the use 
of a Wideband Optical Very Long Base Line Interferometer • 

It uses the very high resolution .intrinsically available with the 
WOVLBI to determine an atmospheric correction in a manner similar to 
the method used in the two-color laser *hiistance measurements ♦ The. 
correction in stellar declination is acconplished by measuring the ap- 
parent declination of a star at two different wavelengths. The dif- 
ference in apparent angle beti^een these tr-jo measurements is then used 
to determine the magnitude of the total refraction correction* The 
accuracy of the declination measurement which results from the use of 
this procedure is by a factor of 30 less than the intrinsic angular 
resolution* However the extremely high angular resolution available 
with a WOVLBI system permits an accuracy of 20 cm or better on the 
polar position* 

In general, when the declination of a star is measured, the 
refraction in the atmosphere results in an offset of the measurement* 

Tor observations made near the zenith, measurements of the local 
pressure and its horizontal gradiant may be used in a quasi-static 
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model of the atmosphere to Tnake corrections to about 0".005. However, 
dynamic effects in the atmosphere (winds)- may cause deviations from 
the quasi— static model which, at the worst type of observing site, 
could he as large as 0".l. 

However, the extremely high resolution (0”. 00008) of the WOVLBI 
makes available a method which does not require the meteorological 
monitoring. To illustrate this method, let us consider the index of 
refraction of the atmosphere as a function of wavelength, for a par- 
ticular pressure and temperature. This may be represented by 

n(X) -1 = 0.000288 |l “(t/ ] 

This expression may be used to obtain the difference in the index of 
refraction for two separate wavelengths* This may be approximated 
by (where 

(^Y) 

- !2l(„(X^).l)+(„a^)-l)] '^X^/ \\) J 

The choice of the critical short wavelength, X , (which essentially 

a 

determines the degree to which the resolution is degraded in obtain- 
ing the refraction correction) is limited by the absorption of the^ 

o 

atmosphere. Thus we choose 3500 A as the center of the short wave- 
length passband. For the long V 7 avelength passband, the limit is 
set by the decreasing sensitivity of the photocathode of the photo- 
multiplier. However, correction factor is rather insensitive to 

© 

this choice j and we will presume a passband centered at 6500 A, so we 


have 



3 


A ^ A(3500 A, 6500 A) - 0.0324 

\ a / 

The above expression describes the properties of the atmosphere. 
Xet ns now relate this to the WOVLBI measurements and their correction. 
Basically, a \XSI-type men sur emeu*: (either radio or optical) measures 
the range to the object of interest from two different reference points 
on earth. While these measurements are normally converted to an angle 
for discussion, the analysis of the atmospheric effects may be more 
effectively carried out in terms of the range. 

Between each reference point and the star, the light passes 
through a part of the atmosphere. Thus for the path to each reference 
point, there is an atmospheric contribution to the range due to the 
fact that the index of refraction of the atmosphere is net equal to 
linity. This deviation in the path length, and for the *tA 7 o 

reference points is given by 

«0 
* 

^J. " d£ 

o 

vjhere n depends upon £ due to the chanses in density, temperature and 
composition alo'ng the path from reference point 1 to the star! If we 
a'ssuma that the variation of composition along the path has a negligihle 
effect (the prii«iry component, water vapor, indeed has a very small 
effect) then the only change in (n - 1) is a multiplicative constant 
and vre may separate the density and vj^velength dependence in the form 
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{n(X) -1} dl 


~ iriiX) -1) 




di = {n(£) -1} R, 


vhere is defined as 




f PiU) 


d5. 


The zero order fringe is deviated from its proper position 
(the position in the absence of the atmosphere). Thus for one par- 
ticular wavelength, the offset, A£(X),^ue to the atmosphere is 

A£(a) = " fn(X) -1] (R^ - R^) . 

Now considering the zero fringe position for tV70 different wavelengths, 
we have, for their separation 

6A£ = A£(X ) ^ A£(X.) = {tn(X ) -1] - [n(X, ) -1]} (R^ ~ R„) 

a D a o , jl j. 

Knowing n(X), we may solve for (Rj^ — R^) a.nd define 




A£(X) = [n(X)-l] [6A£]/[(n(X^)-D - (n(X )-!)] = 6A£/A 

a . D 


we have 
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The quantity 6AZ is the observationally measurable quantity, the 
linear difference between the positions of zero -phase difference. 

One then multiplies by 31 to obtain the total correction required. 

In the same manner, the uncertainty in determining Shi or the un- 
certainty in measuring the position of the zero fringe is multiplied 
by 31 so the accuracy is less than the intrinsic instrumental 
■resolution. 


Douglas G. Currie 
27 December 1971 
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